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In this paper we demonstrate how the recently reported exactly and quasi-exactly solvable models 
admitting quasinormal modes can be constructed and classified very simply and directly by the 
newly proposed prepotential approach. These new models were previously obtained within the 
Lie-algebraic approach. Unlike the Lie-algebraic approach, the prepotential approach does not 
require any knowledge of the underlying symmetry of the system. It treats both quasi-exact and 
exact solvabilities on the same footing, and gives the potential as well as the eigenfunctions and 
eigenvalues simultaneously. We also present three new models with quasinormal modes: a new 
exactly solvable Morse-like model, and two new quasi-exactly solvable models of the Scarf II and 
generalized Poschl- Teller types. 



I. INTRODUCTION 

Quasinormal modes (QNM) has attracted great interest in recent years They arise as waves emitted by a 
perturbed neutron star or black hole that are outgoing to spatial infinity and the event horizon. Generally, the wave 
function of QNMs has discrete complex frequency, whose imaginary part leads to a damping behavior. QNM carry 
information of black holes and neutron stars, and thus are of importance to gravitational-wave astronomy. In fact, 
these oscillations, produced mainly during the formation phase of the compact stellar objects, can be strong enough 
to be detected by several large gravitational wave detectors under construction. 

As black hole potentials are generally too complicated to allow analytic treatment, so in order to understand the 
origin of the discrete imaginary frequencies, one can try approximating the top region of the black hole potential by 
some inverted potentials which are solvable. This has been done by using the inverted harmonic oscillator [2], and 
the Poschl- Teller potential ||. 

Recently, in Q we have extended the number of exactly solvable models that admit QNMs. We take QNMs to 
include both decaying and growing modes with complex energies. Furthermore, we have provided the first model 
with QNM that is quasi-exactly solvable (QES). A system is called QES if a part of its spectrum, but not the 
whole spectrum, can be determined analytically [oT-[l^|. Our approach inJ3|was to study solutions of QNM based 
on the s/(2)-Lie- algebraic approach to one-dimensional QES theory 0, 0, H/Q- We demonstrated that, by suitably 
complexifying some parameters of the generators of the sl(2) algebra while keeping the Hamiltonian Hermitian, we 
could indeed obtain potentials admitting exact or quasi-exact QNMs. These models were later re-studied numerically 
by the asymptotic iteration method in [13j . 

In this paper we would like to show that exactly solvable and QES models with QNMs can be constructed much 
more simply without resorting to the machinery of Lie-algebra. This is achieved through the prepotential approach 
proposed recently p^4l7j . This is a simple constructive approach, based on the so-called prepotential @, EH EH > which 
can give the potential as well as the eigenfunctions and eigenvalues simultaneously. The novel feature of the approach 
is that both exact and quasi-exact solvabilities can be solely classified by two integers, the degrees of two polynomials 
which determine the change of variables and the zero-th order prepotential. Hence this approach treats both quasi- 
exact and exact solvabilities on the same footing, and it provides a simple way to determine the required change of 
variables, say x, to a new one z = z{x). All the well-known exactly solvable models classified in supersymmetric 
quantum mechanics (SUSYQM) [2(|, the QES models discussed in [a. Irjl. lll-tToj . and some new QES ones (also for 
non- Hermitian Hamiltonians), can be generated by appropriately choosing the two polynomials. Our approach, unlike 
the Lie-algebraic approach, does not require any knowledge of the underlying symmetry of the system. Furthermore, 
our approach can generate the Coulomb, Eckart, Rosen-Morse type I and II models [Iff which are not covered by 
the standard Lie-algebraic program yl Hj. Compared with SUSYQM, our approach has the advantage that we do 
not have to assume the sufficient condition for integrability needed in SUSYQM, namely, shape invariance. In fact, 
shape invariance comes out automatically from this approach (l7j ]. What is more, the transformation of the original 
variable z{x) is determined within the prepotential approach, whereas in SUSYQM this has to be taken as given from 
the known solutions of the respective models before one could solve the shape invariance condition. 

We shall adopt the prepotential approach here to generate and classify all one-dimensional exactly solvable and 
QES models with QNMs based on sinusoidal coordinates, namely, those coordinates z(x) whose derivatives (with 
respect to x) squared are at most quadratic in z. Our strategy is to suitably complexify some or all of the parameters 
in the prepotential while keeping the resulted potentials real. We find that all the systems reported in Q can be very 
easily constructed. During the course of investigation, we also realize that there are three new QNM models which 
were missed in Q. We thus take this opportunity to report on them. 
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The plan of the paper is as follows. In Sect. II we briefly review the essence of the prepotential approach. Sect. Ill 
to V then discuss the prepotential construction of QNM models of the Scarf II, the Morse, and the generalized Poschl- 
Teller type, respectively. These three types of potentials have some interesting features, and serve as good examples 
to illustrate the procedure. Furthermore, there is one new QNM model in each of these three types that we would 
like to report. Sect. VI summarizes the paper. In Appendices A to C, we list all other cases for easy reference. 



II. PREPOTENTIAL APPROACH 

The main ideas of the prepotential approach [Til Il5j | can be summarized as follows (we adopt the unit system in 
which h and the mass m of the particle are such that h = 2m = 1). Consider a wave function 4>n(x) (N: non-negative 
integer) which is defined as 

<f> N (x) = e- w °^p N {z), (1) 

with 

PN{z) = \TlL^-^ n>o. (2) 

Here z = z(x) is some real function of the basic variable x, Wq(x) is a regular function of z(x), and z^s are the 
roots of pn(z). The variable x is defined on the full line, half-line, or finite interval, as dictated by the choice of z(x). 
The function pn{z) is a polynomial in an (N + l)-dimensional Hilbert space with the basis (1, z, z 2 , . . . , z ). Wo(x) 
defines the ground state wave function. 
The wave function (f>N can be recast as 

4> N = cxp (-W N (x, {z k })) , (3) 

with Wn given by 

N 

W N (x,{z k }) = W (x)-^2bi\z{x)-z k \. (4) 
fc=i 

Operating on <f>N by the operator —d jdx 2 results in a Schrddinger equation Hn4>n = 0, where 

Bn = ~ + V n , (5) 

V N = W$-W£. (6) 

Here and below the prime represents derivative with respect to x. Since the potential Vn is determined by Wn, we 
thus call Wn the iVth order prepotential. From Eq. one finds that Vjv has the form — V + AV^- 

V = W' 2 -W(,', 

^'-T)§^ + g ( 2 -4-,) - (7) 

Thus the form of Vjv , and consequently its solvability, are determined by the choice of Wo (x) and z' 2 (or equivalently 
by z" — {dz 12 /dz)/2). Let Wqz' = P m (z) and z' 2 — Q n (z) be two polynomials of degree m and n in z, respectively. 
The variables x and z are related by 

and the prepotential Wq(x) is determined as 



Wo(x)=[ I dz^r) . («)) 



z(x) 



We assume ([8]) is invertible to give z = z(x). Eqs. ([8} and (|9]) define the change of variables z{x) and the corresponding 
prepotential Wq(x). Thus, P m ( z ) and Q n {z) determine the quantum system. Of course, for bound state problems 
the choice of P m and Q n must ensure normalizability of cf>o — exp(— Wo)- 

Now depending on the degrees of the polynomials P m and Q n , we have the following situations [Til. [l5j: 
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(i) if max{m, n— 1} < 1, then in Vn(x) the parameter N and the roots Zfc's will only appear as an additive constant 
and not in any term involving powers of z. Such system is then exactly solvable; 

(ii) if max{r7i, n — 1} = 2, then N may appear in the first power term in z, but Zfe's only in an additive term. If N 
does appear before the z-term, then the system belongs to the so-called type 1 QES system defined in i.e., 
for each N > 0, Vjy admits N + 1 solvable states with the eigenvalues being given by the N + 1 sets of roots 
Zfc's. This is the main type of QES systems considered in the literature; 

(hi) if max{m, n — 1} > 3 , then not only N but also z^'s may appear in terms involving powers of z. If Zk's do 
appear before any z-dependent term, then for each N > 0, there are N + 1 different potentials Vn, differing 
in several parameters in terms involving powers of z, have the same eigenvalue (when the additive constant, or 
the zero point, is appropriately adjusted). When Zfc's appear only in the first power term in z, such systems are 
called type 2 QES systems in [y]. We see that QES models of higher types are possible. 

This gives a very simple algebraic classification of exact and quasi-exact solvabilities. Previously exact and QES 
systems were treated separately. 

In the rest of this paper we shall consider only cases with m, n < 2. Coordinates with n < 2 are called the 
sinusoidal coordinates. Let P 2 (z) = A 2 z 2 + A\z + A and Q 2 (z) — az 2 + j3z + 7. Hence the solvability of the system 
is determined solely by A 2 : exactly solvable if A 2 = 0, or (type 1) QES otherwise. The potential Vn takes the form 



V N = W? -W{ + aN 2 -2A X N -2A 2 Nz-2A 2 ^z k -2^^—\p 2 {z k )-%Zk-j -^T^^'' 

fe=i fe=i z Zh I 



2 4 <— ' z fe - zi 



(10) 



Demanding the residues at Zfe's vanish gives the Bethe ansatz equations satisfied by the roots Zfc's: 

a p ^ Qa(zfc) 



p2{Zk) ~ -zZk ~ 7 - V =0, fc = l,2,...,JV, 

/ ZL ^ — ^ Vi. — r r 
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(11) 



A 2 z2 + ( At - -) zfe + A - !j - fc ^ = 0. 



4 ^fc Zfc ~ Zi 



(12) 



Using 



we arrive at the potential 



^(z) = -^=4=, <(z) = z' 



dz 



(13) 



W(x) = 



, , lift 1 In <*Q: 



A' 



- 2Ai7V - aN 2 + 2A 2 Nz + 2A 2 ^ 



Zfe 



fc=l 



(A 2 z 2 + A lZ + A 
az 2 + fiz + 7 



- 2(7V + l)Ajz + - (A 2 z 2 + A x z + A 



2az + (3 
az 2 + j3z + 7 



iV 



(27V + l)A x - aN 2 + 2A 2 ^ z k 

k=l . 



(14) 



and the wave function 



iP N ~ e- M/ > N (z) 



Pat(z). 



(15) 



1 We take this opportunity to correct a typographic error in \H 1 15l | , where the condition for m and n for this case was erroneously written 
as min{m, n — 1} > 3. 
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Eq. (fl~4"]) gives the most general form of potential, based on sinusoidal coordinates, that cover both the exactly and 
quasi-exactly solvable systems. 

It turns out that there are only three inequivalent canonical forms of the sinusoidal coordinates [l7| , namely, (i) 
z' 2 = 7 f 0, (ii) z' 2 = f3z (P > 0), and (hi) z' 2 = a(z 2 + S) (S = 0,±1 for a > 0, and S = -1 if a < 0). Case (i) 
and (ii) correspond to one- and three-dimensional oscillator-type potentials, respectively. In case (iii), for a > 0, the 
potentials are of the Scarf II (S — 1), Morse (<5 = 0) and generalized Poschl- Teller (S = —1) types, while for a < 
and S = — 1, the potentials generated belong to the Scarf I type. 

For clarity of presentation, in the main text we shall illustrate the prepotential construction of QNM models only 
for the Scarf II, Morse, and generalized Poschl- Teller type potentials. Other cases are summarized in the Appendices. 
As mentioned in Sect. I, we choose to discuss these three cases because they have some interesting features that serve 
as good examples to illustrate the procedure, and because there are new QNM models in these types not realized in 
1- 



III. SCARF II: z- 



a{z 2 + 1) 



Consider first the case z' 2 = a(z 2 + 1) (a > 0), which is solved to give z(x) = ± sinh(y / aa;). For definiteness we 
shall take z(x) — sinh(y / aa;). The case corresponding to the negative sign is simply the mirror image of the present 
case (i.e., by taking x — > —x). The same applies to the other sinusoidal coordinates discussed in the rest of the paper. 

Putting p — and 7 — a in (|14p . we get the potential (—00 < x < 00) 



V N = ^(z 2 + l) 



+A 2 



2A 1 
a 



(A - A 2 f 



27V - 1 z 





»)] 


(-- 







+ (Ao-A 2 )( — + l Jz2 + i 



z 2 + l 

|2 



A 9 A 

2A ± N - aN 2 i -(A a - A 2 ) + 2A 2 V 

a a f— ' 



N 

,Zk 
fc=l J 



(16) 



In terms of x, it is 



A 2 

Vn(x) = — cosh 2 (i/ax) 



(A - A 2 f 



sech 2 (-\/ax) 



' 2A \ /2A 
+ (A - A 2 ) ( — + 1 J tanh(V«a;)sech(v / ^x) + A 2 f - 2N - 1 ) sinh(V«x) (17) 



A2 OA N 

2A X N - aN 2 - -± 2 -(A - A 2 ) + 2A 2 V 

a a f— ' 



fe=i 



The prepotential Wq, obtained from ©, is 

Ax 



A 2 Ax , 9 . A a - A 2 , 

Wo(x) = — z + -i In (z 2 + 1) + — -tan- 1 ^ 

a 2a a 

= — - sinh(v / ax) H lncosh(v / ax) H tan -1 sxxxh(^/ax). 

a a a 



(18) 



Hence the ground state wave function ipQ ~ e is 



ipo ~ e" 



(cosh(y/ax)) 



li A 2 -An 



tan 1 sinh(^/o(x) 



(19) 



A. A 2 = 



Let us first discuss the exactly solvable case, with A 2 = 0. The potential becomes 



V N = 



'^-Ax 






a 





sech 2 (v / ax) + Aq ( h 1 ) tanh( v / ax)sech(v / ax) 



2AxN~aN 



A 2 
2 A x 



(20) 
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Suppose all Aj's are real, the potential Vn is just the Scarf II potential given in J20J. Recall that in our approach we 
have Hn4>n — 0, hence the eigenvalue of Hn with potential Vn is always zero. If we define the Scarf II potential 
V(x) by the first two N- independent terms, then Vn(x) — V(x) — En, where En = 2AiN — aN 2 — A\/a are the 
eigenvalues of V(x). The corresponding eigenfunctions are given by ([T]), with </>o given by (|19|). and pn( z) defined by 
the roots Zf-'s of the corresponding Bethe ansatz equation (|12|). which gives the Jacobi polynomials [22| . 
If we take instead 



d 2A 1 

Aq = -i», 

2 a 



then the potential is 



V N (x) 



1 

4a 



(« 2 



1 = -*-, 



cd 



(21) 



d 2 ) sech (-y/aa;) tanh(-y/ax)sech(v / o7a;) 

2a 



{ N + - 

4a V 2 



(22) 



This is the exactly solvable case 1 hyperbolic QNM system discussed in 0]. The special case where d = has been 
employed in Q to study black hole's QNMs. If we take the potential to be defined by the first two terms, then the 
energies and ground state wave function are 



and 



E N = --[N + ~„ 



N+ 2 



tp {x) ~ (cosh v ^x) (lc+Q)/2a exp(idtan- 1 (sinhVo7x)/2a). 



(23) 



(24) 



Note that En is independent of d, which is a general feature of the Scarf-type potentials. Also, the imaginary part is 
proportional to N + 1/2, which is characteristic of black hole QNMs. 



B. A 2 £ 

If A2 0, then reality of the first term of Vjv in (|17p requires that A2 be real, or purely imaginary. 

If A2 is real, then it is clear from the first term of (|19p that ipo, which governs the asymptotic behaviors of ipN, is 
not normalizable on the whole line. Hence there is no QES model with real energies in this case. 

Suppose A 2 is purely imaginary, say A 2 — ic ^ with real constant c. Then the fourth term of Vjv in ([TTf can be 
real provided that 

2A 

— - — 2N — 1 = id, d: real. (25) 
a 

If d 0, then reality of the third term of Vat demands that Aq — A 2 = ±[2(JV + 1) — id]a. But then, as can be easily 
checked, with these values of Ai's the second term of Vat cannot be real, unless d = 0. 

So we are left with the choice A 2 = ic 7^ and 2A%/a — 2N —1 = 0. This implies that Aq — A 2 must be real for 
Vn real. Let Aq — A 2 = act with real a. The potential of this system assumes the form 



Vn{x) = cosh 2 (i/ax) + a 



"" ~ ' N + v ( N + 1) sech2 (^" x ) 



2aa (N + 1) tanh(A/ax)sech(v / as) 




/V 



2icJ2 z k)- (26) 
fc=i / 

The ground state wave function is 

^ Q ^ e -if smh(v^x) ( C osh( v ^x)) _(A ' + ^ ) e -«tan- 1 si„h( v ^x)_ ^ 



For a ^ 0, the potential (|26| is a new QES model with QNMs which has been overlooked in our previous study 
based on the Lie-algebraic theory of 

The case a — leads to a totally different model with real QES energies, which was first discussed in [2l|. We 
briefly discuss it in the next subsection. 
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C. Singular potential with A2 7^ 

For a — 0, i.e., A 2 = Aq = ic ^ 0, the potential (|26|) becomes 



V N (x) = -^-cosh 2 ( % /ax) - a (^N + ^ ^iV+ sech 2 ( % /ax) - ^-j + 2ic^ 



z k • (28) 



This is a singular potential unbounded from below. Yet it exhibits very peculiar features, such as the existence of 
QES bound states with real energies, and QES total transmission modes. We refer the reader to 21] for a detailed 
discussion of this potential. 

At first sight it may seem strange that the system has QES real energies, owing to the term 2ic^2 k z k . We prove 
below that this term is indeed real, as the sum ^ fc z k is purely imaginary. 
The Bethe ansatz equations (fT2|) in this case are 



7+1 

iczl + aNz k + ic - a V " k = 0, k = 1, 2, . . . , N. (29) 



Multiplying (l2l?f by —1 and taking its complex conjugate, we have 

ic(-z k f + aN(-z k ) + ic- «V - ( J** )2 , + -x = 0, k = 1, 2, . . . , N. (30) 

Here f represent the complex conjugate of z. Comparing (f2Uf and (f3"0)) , we conclude that if z^'s are the solutions of 
(|29p . then so are their negative complex conjugates — Z/-'s. This implies that the sum z k is purely imaginary, and 
hence the QES energies are real. 



IV. MORSE: z' 2 = az 2 (a > 0) 

In this case the change of variables is given by z(x) = exp(±-y/ax). Again we shall only consider the positive case, 
i.e., we take z{x) = ejip(y/ax). The potential is (—00 < x < 00) 

V N{ x) = ^W^-ov-iW^f*^ 

m V m I \ m 1 7 m 7 Z \ mm < ^ 



a J \ a j z a z* \ a a 

+A 2 (^-2N-l)e^ x + A (^ + l) e-^ x + 



fc=l 



( 2A X N - a N 2 -^- + 2A 2 V z k ) . (31) 



The prepotential Wq is 



a a 

fe=i 



W = - ( A 2 z + A 1 \n z- — 
a V z 



^e^ + ^l^-^e-^x (32) 
a a a 



Hence the ground state wave function is 



A. A 2 = 

For A 2 = 0, the potential reads 

V N (x) = A + l) e-^ x + ^ e ~ 2 ^ x - ( 2A X N - a N 2 -^-). (34) 

V a ) a \ a 
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To get a real potential defined by the first two terms of Vn, Aq and A\ must be both real, or both complex. 

If A is real, then the first term of Vn requires that Ai be real. The resulted model is the exactly solvable Morse 
potential. 

If Aq is complex, then the second term of Vn demands that Aq be purely imaginary, i.e., Aq = ic with c real. The 
first term of Vn then requires that 2Ai/a + 1 = — id is also purely imaginary. The potential is 



Vn(x) = cde 



e~ 2 ^ x - a 



d 2 



1 



N 2 - N - id [N + 



1 



with ground state 



if e-^ + i(id+l)v^2 



(35) 



(36) 



This gives a new exactly solvable QNM model. This system was overlooked in [J] . 

We note here that this system, or rather its mirror image, can be obtained by a different complexification of Aq, Ai 
and A 2 . We shall discuss this in the next subsection. 



B. Ai 



If A 2 ^ 0, then A 2 has to be real or purely imaginary. 

If A 2 is real, then all Ai's must be real. In this case, for A 2 > and Aq < 0, the potential Vn in (|3"Tj) defines a 
QES system. 

For A 2 purely imaginary, Ai/a — (N + 1/2) must either be zero, or purely imaginary. If A\ja — (N + 1/2) = 0, we 
can have 



An 



b . 
-1-, Ai 
2' 



N+-)r,. L 



This leads to a potential 



V N (x) 



4q 



e 2VZx _ (iv + l)de-^ x 



4a 



6, d : real. 



a . bd 

4 ~ *2a 



N 



ib ^2 z k 



(37) 



(38) 



fe=i 



The system defined by the first three terms of Vjv is a QES system with complex energies given by the last term in 
bracket. This is indeed the very first QES QNM model presented in 

On the other hand, if A\j a — (N + 1/2) is purely imaginary, then one can choose 



An 



A x d 
-ic, — = i- 



a 2 

with c and d real constants. The resulted potential is 



V N (x) = cde^ x - — e 2 ^ x 
a 



N ■ 



An = 0, 



1 d „ c v ^ 
i 2i- > 



2 A . 



(39) 



(40) 



In this case though A 2 ^ 0, the parameter iV and the roots ZkS do not appear in any re-dependent term, and hence 
the system is exactly solvable. 

This system is indeed the mirror image of the model described by (|35[) , although the functional forms of the two 
potentials look very differently To show this, we need to demonstrate that the potential, the eigenvalues and the 
eigenfunctions of one system are mapped into those of the other system under parity transformation x — > —x. 

The wave function of the present system is 



N 



N 



(id+2N+l)^/ax 



Y[ - 



k=l 



where the roots z/c's satisfy the BAE (from <jX2j) ) 



/ d \ N 2 

a U-+JV z fc -aV^=0, fc = l,2,...,iV 

\ / / * * 7l 77 



2#fc 



Zk - zi 



(41) 



(42) 
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We stress here that z is z — e^° x as before. 

Under parity transformation x — > —x, we have z — > z~ x , Zk — > ZjT 1 . Eqs. (|40)) . (|41|) and (|42|) are mapped, 
respectively, into 



and 



- ic— 7? + a [ i— + N 



cde~^ x e" 2 ^ - a 

a 

ijL e -^*+i(id+2N+l)^/ax 



1 .d „ c \ 1 
--i 2z-> — 

9 /-^ y,. 



n 



i i 

z z k 



Zk 



Multiplying (|4"5)) by — z 2 , and using the identity 



1 N 



Z^fc 



Zk{zi - z k ) 



0, fc = 1,2,..., N. 



(43) 
(44) 

(45) 



N 

E 



JV 



N-1-J2 



Zk 



l^k 



Zk - zi 



(46) 



we can transform (f4"5"]l to 



AT 



a \i- + l zt+ac-aV = 0, fc = 1, 2, . . . , N. 



Zk - zi 



(47) 



This is just the BAE for the roots z^s of the eigenfunctions of the system defined by the potential (|35|). Hence the 
BAE of this system is mapped into the BAE of the system defined by (|35[) under parity. 
Now the eigenfunctions (|44"1) can be rewritten as 



N 



N N 

e'^^+jtW+i)^ (e^™ x ) z~ N JJ ( 



z- Zk), 



fe=i 



N 

n< 

fc=i 



z - z k ) , 



(48) 



where we have used the fact that (e^ x ) N z~ N = 1. Eq. (|48l) is just the eigenfunction of the potential (|3"5|) . Thus, 
together with the result of the last paragraph on BAE, one sees that under parity the eigenfunctions of the system 
given by (14*01) are mapped into those of given by (|"">5j) . 

Finally we show that the last term of the potential (f4T"j) is equal to the last term of ([55]) . This proves the equality 
of the eigenvalues of both systems. Dividing the BAE (|47p by Zk, summing over all k and using the identity 



N 



EE 



Zk 



1 



Zk - zi 



N(N-1), 



(49) 



we get 



N 1 
„ c \ -\ 1 
2«c- > - 



fe=l K 



So the last term of the potential (|40[) becomes 



a „ c \ - 

- z 2i- > 



-zdA — AT — -/V. 



-N 2 -N-id(N+- 



(50) 



(51) 



This is equal to the last term of ([55)1 . Hence under parity potential (14TJ1) is mapped into (f3"5| . 
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V. GENERALIZED POSCHL-TELLER: z' 2 = a(z 2 1) 

In this case z(x) — cosh(^fax) . The potential is (0 < x < oo) 

'(A + A 2 ) 2 



V N (x) 



A 2 

— - sinh (\Tolx) 



At 





»)] 


(-- 


V a 





cosech (\/ax) 



'2A \ /2A 
+ (A + A 2 ) ( + 1 J coth( % /aa;)cosech(\/ax) + A 2 ( — - 2N - 1 ) cosh( % /ax) 



(52) 



9 -4? 2A 2 

2A X N -aN 2 

a a 



N 



(A„+ A 2 ) + 2A 2 ^2z k 

k=l 



5 s - e A2 + A ° coth- 1 cosh( x /Sa:) 



The ground state wave function ipo ~ exp(— Wo) is 

?Ao ~ e"^^^ (sinh(V^a;)) 
~ e -4?™*(^) ( sin h(V^a;)) 
Since Vjv — >■ oo as x — > 0, we must have the boundary condition ipN — > as n- 



tanh(-^x) 



(53) 



For A 2 = 0, the potential is 
V N 



r a 2 








-)] 


a 







A. ,4 2 = 



2A 

cosech 2 (Vox) + A ( h 1 ) coth(v / ax)cosech(v / ax) 



A 2 ' 

2AiiV-aiV 2 - 

a 



(54) 



If all Aj's are real, this gives the generalized Poschl- Teller potential in [20j. If we take the Poschl- Teller potential 
V(x) to consist of the first two N- independent terms, then the eigen-energies are En — 2AiN — aN 2 — A\/a. These 
eigenvalues are the same as those of the Scarf II potential. Our approach makes it very easy to see why the eigenvalues 
are the same for the two apparently different systems. 
We can also take 



A 



d 2Ai 

-V — + 1 

2 a 



(55) 



which result in the potential 

V N (x) = 



4a 



(a 2 + c 2 + d 2 ) cosech (Jax 



cd 
2a" 



coth(y / ax)cosech(v / aa;) 



c 

4a 



N+l) a-zc(A + I 



(56) 



This is the exactly solvable case 2 hyperbolic QNM system discussed in [4J. 

B. A 2 / 

As with the Scarf II case, if A 2 ^ 0, then reality of the first term of Vjsr in (|52|) requires that A 2 be real, or purely 
imaginary. 

If A 2 is real, then Ao and A\ have to be real. The boundary condition that ipN — > as x — > oo is met if A 2 > 0. 
As x — > 0, we have 



ipo — x 



A 2 +A 1 +A 



(57) 
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Thus tpo — > as x — ¥ is guaranteed if 



A 2 + Ax + A < 0. 



(58) 



With AiS satisfying (1581) and A 2 > 0, the potential Vat gives a QES system with real energies. 

If A 2 is imaginary, say A 2 — ic ^ with real constant c, then the fourth term of Vpj in (|52[) can be real if 
2 A x / a - 2N - 1 = or id ^ 0. As in the Scarf II case, if 2A x a - 2N - 1 = id ^ 0, then reality of the third term of 
Vn demands that Aq + A 2 = ±[2(N + 1) — id]a. But then the second term of Vjv cannot be real, unless d = 0. 

So we must have 2Ai/a — 2N —1 = 0. In this case one has A 2 + Aq = —aa ^ with real a. The potential Vjv is 



V N (x) 



■ sinh 2 (\/aa 



N 



N ■ 



cosech (vax 



-2aa (N + 1) coth(- v /aa;)cosech(v / ax) 

N \ 

2ica + 2ic Zk J , 

k=l J 



(59) 



o 
4 



and the ground state wave function 
To satisfy the boundary condition 



■ coshf-yaa:) 



sin 



h(v / ax)) 



-(w+l) 



tanh(^^a;) 



as x — > 0, we have from (|58|) the condition 
a > iV+ -. 



(60) 



(61) 



With these Ai's, we have a new QES system with QNMs. 

Unlike the Scarf II case, however, when a = 0, we do not have a QES model with real energy, since in this case tpo 
is not normalizable at x — 0. 



VI. SUMMARY 



Exactly solvable models admitting QNM maybe useful in providing insights to QNMs emitted from more complicated 
systems such as black holes. In 0] we have found some new exactly solvable QNM models, and a new QES QNM 
model within the sl(2) Lie-algebraic approach to QES theory. 

In this paper we have demonstrated how exactly solvable and QES models admitting quasinormal modes can 
be constructed and classified very simply and directly by the prepotential approach. This approach, unlike the 
Lie-algebraic approach, does not require any knowledge of the underlying symmetry of the system. It treats both 
quasi-exact and exact solvabilities on the same footing, and gives the potential as well as the eigenfunctions and 
eigenvalues simultaneously. We also present a new exactly solvable Morse-like model with quasinormal modes, and 
two new quasi-exactly solvable models with quasinormal modes of the Scarf II and generalized Poschl- Teller types. 



Acknowledgments 



This work is supported in part by the National Science Council (NSC) of the Republic of China under Grant NSC 
96-2112-M-032-007-MY3 and NSC-99-2112-M-032-002-MY3. 



Appendix A: Shifted oscillator: z' 2 — 7 > 



The corresponding transformation z(x) is z(x) — sF/x + constant. By an appropriate translation one can always 
set the constant to zero. Hence we shall take z(x) — y/jx without loss of generality. 
The potential is (—00 < x < 00) 



V N = A 2 2l x A + 2A 2 A lv ^x 3 + {A\ + 2A 2 Ao) x 2 + 2 



AxAo 



A 2 {N+ 1) 



1 A 2 A 
2Ai(N+-)-^-+2A 2 ^2z k 



fc=i 



(Al) 
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The ground state wave function fa ~ exp(— Wo) is 

fa ~ e~^^ x (A2) 

We want Vjv to be real. If A 2 ^ 0, then the term A 2 (N + 1) in the fourth term implies A 2 be real. However, this 
implies the wave function faj, whose asymptotic behavior is governed by fa, is not normalizable on the whole line. 
Hence there is no QES models, with or without QNMs. 

For A2 = 0, the potential is 

V N = (a iX + ^ -2A 1 {N + \). (A3) 

So A\ and A$ must both be real, this is just the potential of the shifted oscillator. 
If we take 

^i = -i~, ^o = 0, (A4) 



the potential becomes 



4 

This is the case 5 exactly solvable QNM model discussed in [4| 



V N = -\c 2 x 2 + ic[N + . (A5) 



Appendix B: Radial oscillator: z' 2 = f3z, (/3 > 0) 



For simplicity we take z(x) = jx 2 . The potential is (0 < x < 00 

\ 2 

64 2 



Vn 



1 ,A 2 2 ^+ l -A 2 A l px i + \ 



A\ + 2A 2 A Q - -A 2 p (AN + 3) 



4A (Ap 1\ 1 

P V P V x 2 



M 
2 



AN + 1 



AA 



N 

2A 2 ^2 z k 



fe=i 



(Bl) 



The ground state fa ~ exp(— Wq) is 



fa ^ x p e ± AlX 32 A2fix 



(B2) 

It can be checked that if A 2 ^ 0, then all Ai's have to be real. Hence there is no QES model with QNMs in this case. 
For A 2 > and Aq < 0, the wave functions ipN are normalizable on the positive half- line. The system is then a QES 
model with real energies. Together with the discussion in Appendix A, one concludes that one-dimensional QES models 
start with degree six, i.e., the sextic oscillator (23|. In fact, the potential (|B1I) with A 2 = 2a > 0, A\ = 2b, A^ = 
and f3 = 4, namely, 

+ 2abx A + [b 2 - (AN + 3)o] x 2 



V N 



a 2 x 6 



(AN + 1)6 + 4a^ 



Zk 



(B3) 



is the very first QES model discussed in the literature 
If A 2 = 0, then (iBil becomes 



AA fA Q , 1\ 1 



/3 V '27 a; 2 

For real j4i and Ao, this is just the radial oscillator. 
Suppose we take 



A x 2N 



1 2A 

2 /? 



The potential takes the form 



/9 = 4, Ai = -2ia, and A = -27. 



W = _ a V + + la (AN + 27 + 1) . 



(B4) 

(B5) 
(B6) 



This is the Case 4 model with QNMs presented in 
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Appendix C: Scarf I: z' 2 = a(l - z 2 

In this case z(x) = sin(^/ax). The potential is 

'(A + A 2 ) 2 



V N (x) 



-4 



-An 



'(Vox) + 



2A 1 



A t 



Ay 



sec 2 (y/ax) + (A + A 2 ) 



2Ai 
a 



2N+1) sin(Vax) - 
with the ground state wave function 



A 2A 

2AiN + aN 2 + ^ + —(A a + A 2 ) + 2A 2 V 
a a * — ' 



1 J ta,n(y/ax) sec(yfax) 

N 



tpo ~ e" 



[cos(^/ax)) 



S 1 e - A 2 + A tanh^ 1 eia(t/ax) 



k=l J 



(CI) 



(C2) 



The system is defined only on a finite interval, usually taken to be y/ax E [— n/2, tt/2]. Hence there are no QNMs. 

If A 2 ^ 0, then all A^s have to be real. The third term of (|C2[) then implies that the wave functions ipN are not 
normalizable. So there are no QES systems with real energies. 

For A 2 = 0, the potential becomes 



V N 



A 



^ + A 1 



Ay 



sec (y/ax) + A Q 



2A 1 



1 I t&ii{^/ax) sec(\/ax) 



2A X N + aN 2 + 



A 2 

2 i A l 



(C3) 



For real A\ and A$ : this is Scarf I potential given in [2C 
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